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l" 1 INTRODUCTION 

• I— I ■ 

S^ ' The cosmic microwave background (CMB) is currently one 
$_( ' of the most powerful tools of cosmology. Its study pro- 

. . 1 vides us with essential information about the origin and 
evolution of the Universe. In particular, a key issue is 
whether the CMB temperature fluctuations follow a Gaus- 
sian distribution, as predicted by the standard inflation- 
ary model. A detection of intrinsic non-Gau ssianity in the 
CMB w ould be a hint of new physics (See iBartolo et alJ 
I2OO4I or iTurok fc Spergel|ll990l and references therein) and 
therefore would have far reaching consequences on our 
current knowledge of the Universe. In addition secondary 
anisotropics (such as the Sunyaev-Zeldovich effects, or the 
Rees-Sciama effect), contaminant astrophysical emissions 
(coming from Galactic and extragalactic foregrounds) and 
systematics can also leave non-Gaussian imprints in the ob- 
served microwave sky. Therefore, it is crucial to perform a 
careful study of any possible detected non-Gaussianity in 
order to understand its origin. 

In the last years there has been a large number of 
experiments measuring the CMB ani sotropic s and polariza- 
tion (DASI, 'Halvers on et alJ l2002l: VSA, LGrainge et al." 

l2003l: CBI, Mason et alJ 2003t: AC BAR. Kuoetal. 

120041 : Archeops, IBenoit et al.l 1200311 . In particular, 
the NASA WMAP satellite has provid ed with high 
quality observations of the whole sky JBennett et alJ 

120038^ . Many works have studied the Gaussianity of 
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We study the power of several scalar quantities constructed on the sphere (presented 
in Monteserin et al.) to detect non-Gaussianity on the temperature distribution of 
the cosmic microwave background (CMB). The test has been performed using non- 
Gaussian CMB simulations with injected skewness or kurtosis generated through the 
Edgeworth expansion. We have also taken into account in the analysis the effect of 
anisotropic noise and the presence of a Galactic mask. We find that the best scalars 
to detect an excess of skewness in the simulations are the derivative of the gradient, 
the fractional isotropy, the Laplacian and the shape index. For the kurtosis case, 
the fractional anisotropy, the Laplacian and the determinant are the quantities that 
perform better. 
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the WMAP data, including the use of w a velet tools 
iVielva et al.l 120 04'. 'M ukheriee fc \Yand I^OOl ICruz et al 
J2005| ICruz eta l. 2006a. ICruz et al.^^l2006bl . iMcEwen et al 
[2005 and Cavon, Jin & Trcastcr 2005), h armo ni c techniques 
iKomatsu et al.. .2003-. Chiang et alj IJOol IColes et alJ 



200^^^_^NaselskY jDoroshkevich fc Verkhodanovl 
Maeueiio fc Med eirosI |2004^. a n d rea l space 



1200- 
analy- 



ses JParkI l2004l. lEriksen et all l2004al. iHansen et al. 



worKS on reiai eu ineunou s mcmue 

pr opert ies iJBon^^ EfstathimJ 

l" 1997, Gurzadvan fc TorrS 



| 2004|. iLarson fc Wandelti Eooi iLarson fc Wandelt 
I2OO5I ICabella et alj |2005|), finding in some cases devi- 
ations from Gaussianity whose origin is uncertain. In this 
paper, we are interested in the study of the Gaussianity of 
the CMB using a set of scalar quantities of the temperature 
field, such as the Laplacian, the shape index or the Gaus- 
sian curvature. Previous works on related methods include 
the study of ma> 
11987!, "Barreiro ot al." 

il997i . iBarreiro, Martinez-Gonzalez fc Sanz 2001), the 
fracti o n of lake and h i ll poin ts JDore. Colombi fc Boucheti 
I2OO3I ICabeUa et alJ l2005^ or the skeleton l enght 
jNovikov. Colombi fc Dorj|2006llEriksen et al.ll2004bl) . 

Some previous works have also compared the perfor- 
mance of different statistics to detect non-Gaussianity in 
the CMB usin g a reference s e t of n on-Gaussian simulations. 
For instance, iHobson et alj (|l999|) compared the capabil- 
ity of wavelet techniques and Minkowski functionals as non- 
Gaussianity detectors using simulated maps of the Kaiser- 
Stebbins effect. The comparative performance of the Spher- 
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ical Haar Wavelet and the Spherical Mexican Hat Wavelet 
has been tested on non- Gaussia n simulations generated 
through the Edgcworth expansion ( M artinez-Gonzalez et alJ 
[20023. Aghanim et al. ( 2003) studied the relative perfor- 
mance of wavelet and Fourier based methods by applying 
them to three different sets of non-Gaussian maps that were 
representative of possible non-Gaussian signatures present 
in the CMB. Three different approaches i n pixel , harmonic 
and wavelet space were tested by Cabella et alJ (|2004j on 
simulations of the quadratic potential mo del and on Gaus- 
sian maps contaminated by point sources. Ijin et al.l ll2005r 
have compared theoretically and experimentally the perfor- 
mance of the kurtosis to the Higher Criticism and MAX 
statistics in the wavelet and curvelet spac e s. 

In a previous work, iMonteserin et all ll2005ll (hereafter 
M05) studied the probability distribution of a series of scalar 
quantities assuming that the CMB is a homogeneous and 
isotropic Gaussian field and including also the effect of 
anisotropic noise and the presence of a mask. These quan- 
tities were proposed as promising tools for the analysis of 
the temperature distribution of the CMB. In the present 
paper, we study the performance of these scalars (as well as 
three new ones) to detect non-Gaussianity in the CMB. In 
particular, we have tested the method using non-Gaussian 
simulations with skewn ess and kurtosis injected throug h the 
Edgeworth expansion (Martine z-Gonzalez et alJl2002f) . Al- 
though Edgeworth simulations do not correspond to a real- 
istic model, it is expected that physical motivated models, as 
well as contaminant emissions, will introduce a certain level 
of skewness and/or kurtosis in the CMB distribution. As 
an example, the quad ratic potential m odel would introduce 
skewness in the CMB iKomats u fc Sperpjcl 200 1 ^ ). the cosmic 
string s would generate ku r tosis llHobson. Jones fc Lasenbvl 
I1999I . iBarreiro fc HobsonI l200l[) whereas residu al point 
sources would introduce both IIArgiieso et al.ll2006r) . There- 
fore the use of non-Gaussian simulations constructed 
through the Edgeworth expansion, where a certain level of 
skewness or kurtosis is injected, provides us with an inter- 
esting tool in order to test the type of deviations from Gaus- 
sianity that one would expect in realistic data. 

The article is organised as follows. In §2 we briefly sum- 
marise the work of M05, describing the different considered 
scalars (including three new quantities) and presenting some 
theoretical results for the case of a Gaussian field. Our simu- 
lations are described in §3 whereas §4 presents the results of 
our Gaussianity test. Section §5 summarises our main con- 
clusiones. Finally, we include a series of appendices with the 
description of the three new quantities as well as analyti- 
cal expressions for the cumulative functions of some of the 
considered scalars. 



2 SCALARS ON THE SPHERE 

The temperature anisotropies of the CMB are usually de- 
scribed as a 2-dimensional field T{9, (j)) on the sphere. The 
first and second derivatives of the field encode very in- 
teresting information about the fluctuations. In particu- 
lar, different quantities which are scalars under a change 
of the coordinate system (i.e. regular general transforma- 
tion (a;i,a;2) —> {x'i,x'2) : s'{x'i,x'2) — s{xi,X2)) can be con- 
structed from the first and second covariant derivatives of 



the field, which can be useful to perform Gaussian studies 
of the CMB. Since a field can deviate from Gaussianity in 
an infinite number of ways, it is a non-trivial task to design 
a unique set of scalars that is optimal to detect any possible 
type of non-Gaussianity. Therefore our aim is to study the 
performance of a (relatively large) number of scalars to de- 
tect some generic types of non-Gaussianity (non-zero skew- 
ness and kurtosis). In particular, M05 studied the probabil- 
ity density function of different scalars for a homogeneous 
and isotropic Gaussian field. In this section we briefly de- 
scribe these quantities (for a more detailed description, see 
M05) and introduce three new scalars (curvedness, fractional 
isotropy and fractional anisotropy). 

Using first derivatives we can construct the square of 
the modulus of the gradient g, which provides information 
on how the temperature varies spatially, and is given by: 



g = T'X^ 



(1) 



Other scalars can be constructed using only second deriva- 
tives. These scalars are directly related to the eigenvalues A^ 
of the negative Hessian matrix of the field, [—T-ij]- In par- 
ticular the eigenvalues can be written as a function of the 
covariant second derivatives in the following way: 



Ai = -- 



(T-J - ^{r\Y-2{r'\T- 



'Y'Vj 'T'V' 



M 



M + /(r^ 



2 T'lT'- 



T\T''' 



(2) 



,(3) 



Combining these three quantities, we can construct other 
scalars that enhance different features or properties of the 
analysed temperature field. An interesting scalar is given by 
the trace of the Hessian matrix (the Laplacian) A+: 



A4 



-Al — A2 



(4) 



The distortion A_ and the shear y are related to the differ- 
ence between eigenvalues: 



A- = Al — A2 



y 



-{Xi-\2Y 



(5) 
(6) 



While these scalars describe the asymmetry of the field 
through the difference of the eigenvalues, other scalars ex- 
press a similar information using their ratio. The ellipticity 
e and its bounded construction the shape index l are dimen- 
sionless scalars given by: 



Al — A2 
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(7) 



(8) 



Other scalars, related to the curvature of the field, are the 
determinant d of the Hessian matrix, the Gaussian curvature 
KG and the curvedness c (see appendix El: 



d- 
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(9) 
(10) 

(11) 
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We can also construct new scalars by combining two of 
the previous scalars. Examples of these are the fractional 
isotropy fi and the fractional anisotropy fa, which cor- 
respond, respectively, to the Laplacian and the distortion 
weighted by the curvedness (see appendices [51 and IHl : 

f,^J-^l±M=, (12) 



fa = 



1 Ai - Aa 



(13) 



Finally, we have also considered the derivative of the 
square of the gradient modulus, Dg-. 



Da 






(14) 

Assuming that T{0, (j)) is a homogeneous and isotropic 
Gaussian random field, we can derive analytical or semi- 
analytical expressions for the probability density function 
(pdf) of the previous scalars (see M05 and appendices 1X1 to 
1^ . In this case, the pdf of the different scalars is completely 
determined by the power spectrum, C'l of the original field 
T. This dependence appears in the distribution function of 
the scalars through the moments of the initial field, a^: 



-^ = ^^^ 4. 



[£{i + l)]' 



(15) 



However, as shown in M05, it is possible to remove the de- 
pendence of the pdf of the scalars on the power spectrum of 
the field by constructing new quantities that we have called 
normalised scalars. Therefore, for the Gaussian case, a given 
normalised scalar will have the same distribution function 
independently of the power spectrum of the field. The defi- 
nitions of the normalised scalars - as a function of the ordi- 
nary scalars - as well as their pdf's for a homogeneous and 
isotropic Gaussian field are given in table0(for completeness 
we also include the normalised temperature in the table). In 
addition to its greater simplicity, these new quantities al- 
low us to deal in a straightforward way with the presence of 
anisotropic noise (see below). 

We have to point out that many of the considered nor- 
malised scalars have some level of correlation between each 
other. For a Gaussian field the maximum number of uncor- 
related scalars is three, for scalars constructed from the first 
and second derivatives. In particular, the square of the mod- 
ulus of the gradient is uncorrelated with all the other scalars. 
Therefore a set of three uncorrelated scalars will be formed 
by g and two other scalars, whose choice is not unique. All 
the information about correlations between the normalised 
scalars is given in table 2 of M05 and table IXTl of this work. 

In order to construct the normalised scalars of a field, a 
necessary step is the calculation of the moments (tq, ai and 
(72 given by equation H15|l through the Ce's. For a noiseless 
CMB map this is trivial provided that we know its power 
spectrum. However, in the presence of anisotropic Gaussian 
noise and a mask some further considerations need to be 
made. Given that the noise and the mask introduce discon- 
tinuities in the CMB map, which represents a problem for 
the calculation of the derivatives, we first filter the consid- 
ered map with a Gaussian beam of dispersion ag (which is 
chosen to be equal to the dispersion of the beam of the ob- 
served CMB map) . The anisotropic noise is characterized by 
a different dispersion at each pixel, an (x) (with x the unity 



vector on the sphere in the direction of observation). In order 
to construct the normalised scalars we need to calculate the 
moments of the field taking into account the different noise 
in each pixel. This is done using a pixel dependent power 
spectrum Hi{x) (see M05): 



Hi{x) 



c;e- 



t(t+i)al . 4,Tval{x) 



+ 



N^^ 



-e(e+i)a 



(16) 



p^x is the total number of pixels of the map. He{x), 



where Nj 

for fixed x, corresponds to the power spectrum of a map 
containing the filtered CMB signal plus isotropic noise with 
dispersion a„ (x) and filtered again with a Gaussian beam of 
dispersion ag. Note that, due to this extra smoothing, the 
dispersion of the noise in a given pixel will depend not only 
on the noise level on that position but also on the a„ of 
its neighbouring pixels. However provided that (t„(s) varies 
smoothly, as it is usually the case, this effect is small. 

Therefore, using Hi(x), we can obtain the moments of 
the field ai{x) at each pixel and, from them, construct the 
normalised scalars by introducing these pixel dependent mo- 
ments on their corresponding definitions. The normalised 
scalars constructed in this way for a Gaussian initial field will 
follow the pdf's given in table Q which are independent of 
the underlying power spectrum and the level of anisotropic 
noise. 

We would like to point out that if the dispersion of 
the noise does not vary smoothly along the map, equation 
1161 1 may not be a sufficiently good approximation to obtain 
ai{x). In this case, the contribution of the noise to the mo- 
ments of the field can still be obtained using simulations. 
Note that since signal and noise are independent, their con- 
tributions to ai{x) can be calculated separately and then 
simply added together (quadratically) . In order to obtain 
numerically the contribution of the noise to ai{x), we need 
to generate a large number (1000) of noise simulations, that 
are smoothed with a Gaussian beam of dispersion ag. We 
obtain then the dispersion of the noise ctq in the smoothed 
map at a given pixel, as the dispersion of the values of the 
1000 simulations at that pixel. To obtain cri(af) and a2{x), 
we need to construct two ordinary scalars, the square of 
the modulus of the gradient g and the Laplacian A+, cor- 
responding to the previous smoothed noise simulations. As 
shown in M05, the mean value of g at each pixel is given by 
ai{x) whereas 0-2(3;) corresponds to the dispersion of A+ at 
each position. Therefore, by estimating these quantities at 
each pixel from the ordinary gradient and Laplacian of the 
smoothed noise simulations, we obtain the contribution of 
the noise to the moments of the field. 

By comparing the estimated values of ai{x) using the 
numerical and analytical (approximated) methods, we can 
quantify which is the error introduced in these quantities 
when using the approximated technique. In particular, we 
have calculated ai{x) using both methods for the example 
considered in this work, which corresponds to the level of 
noise expected in the Q-l-V-l-W combined map of the 4-year 
WMAP data (see ^ . We find that the average error^ on the 
calculation of the noise contribution to ai {x) is less than 2 



^ This error is defined as 100 X cr" — cr" /ct", where ct" and <t° 
refer to the ith moment obtained with the numerical and approx- 
imated approaches repectively 
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Normalised scalar Notation 



Definition 



Domain 



pdf 



temperature 



T 



{— oo, oo) 



p(^) = ^-' 



greatest eigenvalue 



Al i(A+ + A_) (Aa.oo) 



p(Ai) 



'1 1 



l + er-/( ^§A] 



lowest eigenvalue 



A2 



i(A+-A_) (-oo.Ai) 



P(A2) 



av^ 



^1 1 



'27r ■ 
3 ' 



1 - erf 



Laplacian 



(—00, 00) 



P(A+) 



V2^ 



determinant 



A1A2 {—00,00) 



p{d) 



_4_ 4d 



1- erf [ V 6d 



d < 
J> 



shear 



(0,00) 



p{y) = Ae-'^y 



distortion 



V^ 



(0,00) 



p(A_) = 2A_e 



ellipticity 



xA 



(—00, 00) 



p(e) = 4|e| (l + 8e2) ^ 



shape index 



— arctan 



(-A) (-1.1) 



p(^') = f 



I cos (-It)! 



^[^+2 



COt^ 5? 



fractional anisotropy 



/a 



v^ 



v^ 



+^i 



(0,1) 



P(/a) 



2/a 
(1-/I)*(l + /|)' 



fractional isotropy 



/. 



1 A1+A2 
72 



V^ 



+ A^ 



(-1,1) 



Pifi) 



(2-/7)5 



gradient 



(0,00) 



p(g) = e-5 



derivative of gradient D„ 



ia?yi 



(—00, 00) 



P(^«) = T^/o 



Gaussian curvature kq 



(1 + 9)^ 



{-00,00) p{kg) 



/: 



G _2e_ 

00 KG 



3Sg 



~G e-*^d2 



KG < 



/; 



2e 
G ''G ■ 



^e V ~Ge4^ [i_er/(V6i)] d2 kg > 



Curvcdness 



c |\/^^ + ^- (°'°°) 



p{e) = ^\c\e-^^' D {V2e) 



Table 1. Definition of the normalised scalars as a function of the ordinary ones. The pdf's correspond to a Gaussian initial field and 
are independent of its power spectrum. Note that the values of ct^ may be pixel-dependent. 



per cent for the three moments, with maximum errors less 
than 15 per cent for every pixel. Moreover, if we consider 
the error in (Ji{x), including the contribution of both the 
signal and the noise, the average error becomes negligible 
(less than 0.1 per cent) whereas the maximum error is less 
than 1 per cent for every pixel and for the three moments 
considered. Therefore, for the considered work, we have used 
the simpler approach given by equation llfcill to estimate the 
moments ai (x) since the approximation works very well. 

Regarding the effect of the mask, the smoothing of the 
masked map (with the pixels of the mask set to zero) reduces 
the discontinuty at its boundary. However, this smoothing 
strongly contaminantes the pixels close to the boundary of 
the mask, and therefore the scalars calculated at those po- 
sitions, which must be removed from the analysis. In order 
to do this, we will consider in our study only those pixels 
outside an extended mask that excludes also those pixels 



significantly contaminated by the smoothing of the original 
mask. 



This extended mask is obtained following the heuristic 
approach proposed by (Erikscn ct al. 2004b). First, the orig- 
inal mask is smoothed with a Gaussian beam of dispersion 
3 times the pixel size. Then, all those pixels above 0.991 
are set to 1 (kept from the analysis) whereas those below 
this value are set to (excluded from the analysis). Note 
that this way of constructing the mask is different and more 
conservative than the one presented in M05. In any case we 
have tested that both methods provide very similar results, 
but we have chosen the one of the present work due to its 
simplicity. 
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3 THE SIMULATIONS 

In order to test the performance of the considered scalars 
to detect non-Gaussianity, we have generated different sets 
of Gaussian and non-Gaussian CMB simulations. In par- 
ticular, we have used the Edgeworth expansion to simu- 
late non-Gaussian CMB maps with a certain level of skew- 
ness or kurtosis. Although the simulations obtained in this 
way do not correspond to a particular physical model, they 
are, nonetheless, a useful tool to mimic some generic de- 
viations of Gaussianity expected in realistic data. Indeed 
many physical motivated models, as well as contaminant 
emissions, will produce a certain level of these higher order 
moments on the CMB observations such as the quadratic 
potential infla tionary model (which introduces skewness, 
iKomat su fcS]3ergcl_20 0I), co s mic strings (which genera tes 
kurtosis. iHobson et al.lll999l. iBarreiro fc Hobsonll200ll) or 
point source residuals (which produce both, lArgiieso et alJ 
I2006D . 

The Edgeworth expansion of a one-point density func- 
tion f{y) can be expressed in terms of the H ermite polyno- 
mials (see e.g. iMartinez-Gonzalez et al.ll2002l and references 
therein). Keeping only the first terms in the correspond- 
ing Hermite polynomials and considering only the skewness 
and kurtosis perturbations, we approximate the distribution 
function of our non-Gaussian simulations by the following 
two equations: 



fs {z) = 



fK\ 



{l + f[^(^^-3)]} 



24 ^ 



6z' 



+ 3) 



(17) 



(18) 



where S and K denote skewness and kurtosis, espectively. 
These distribution functions are not always well defined, be- 
cause they can become negative even for relatively small val- 
ues of 5* or 7^. To avoid this problem we set the function to 
zero when it becomes negative and then renormalise it to 
unit area. Note that for values oi S, K < 1, the zeroes of 
the function always appear in the tails of the distribution, 
so the renormalization value is close to 1. 

Our test CMB simulations have been produce d with 
the aid of the HEALPix package JGorski et al.ll2005r . using 
a resolution of Nside = 256 (whi ch corresponds to a pixel 

of 13.7 arcminutes). Following iMartinez-Gonzalez et all 

||2002), we have generated non-Gaussian Edgeworth sim- 
ulations as follows. First we produce simulations of white 
noise with a distribution given by equations 1171 and I18II 
for different values of S and K (considering also the case 
S = K = G, which corresponds to a Gaussian distribu- 
tion) . These maps are then convolved with a Gaussian beam 
of FWHM=23 arcmin. Finally we renormalize the power 
spectrum of the resulting map to the desired CMB power 
spectrum. In particular, we have used the power spectrum 
iven by the best-fi t model to the 1-year WMAP data 
ISpergel et alJ l2003ll convolved with a Gaussian beam of 
FWHM=33 arcmin. The C^ 's for this model were generated 
using CMBFAST ( Seliak fc Zaldarriaga 1996) . 

Note that as a consequence of the first beam convolu- 
tion and the introduction of correlations in the temperature 
maps, the simulated field does not longer follow the expres- 
sions given by equations 11711 and 1181 and, in particular. 



the levels of the injected S and K have been significantly 
reduced. Nevertheless, the important point is that this pro- 
cess generates non- Gaussian simulations with the desired 
power spectrum and whose final level of S and K can be 
controlled through the values of the original injected skew- 
ness and kurtosis. The upper part of table|21gives the mean 
and standard deviation of the final skewness and kurtosis 
obtained from 1000 simulations for different injected values 
of S and K. Note that the reduction in the level of K is 
more important than the one in S. 

As a first step to test the power of the scalars to de- 
tect non-Gaussianity we have applied our method to these 
three sets of ideal simulations: Gaussian, injected skewness 
(5" = 0.08) and injected kurtosis {K=QA) (sectionlCT . Note 
that we have chosen to inject in the simulations only skew- 
ness or kurtosis to discriminate better the effect of each of 
these higher moments in the Gaussianity test. We remark 
that the particular values of 5* and K have been selected to 
allow for a good comparison between the performance of the 
different scalars (although consistent results were obtained 
using different values of these parameters) . 

However, in a realistic experiment, the data will be con- 
taminated by some level of instrumental noise. Moreover, 
some regions of the sky may be severely contaminated by 
different astrophysical emissions and thus will be excluded 
from the analysis. Therefore, we have also analysed a set 
of realistic simulations that take into account these prob- 
lematics ("section 14.211 . In addition to the CMB signal, our 
realistic simulations contain anisotropic Gaussian noise at 
the level expected in the Q-l-V-l-W combined map of the 
4-year WMAP data and ha s been masked using the KpO 
mask of the WMAP team llBennett et alj|2003d) . As ex- 
plained in the previous section, in order to reduce the dis- 
continuities produce by the noise and the mask, we have 
further smoothed the simulations with a Gaussian beam of 
FWHM=33 arcminutes. Finally, following the procedure ex- 
plained in section l^l we have constructed an extended mask 
that excludes from the Gaussian analysis those pixels con- 
taminated by the original mask due to this extra filtering. 
Note that whereas the KpO mask covers approximately a 
23.9% of the full sky, the extended KpO mask used in the 
Gaussian analysis excludes a total of 31.7% of the pixels. 

For the realistic case, we have considered also three dif- 
ferent sets of simulations: Gaussian, injected skewness and 
injected kurtosis. Obviously, due to the presence of noise 
and mask, it is more difficult to distinguish between dif- 
ferent models. Therefore, in order to obtain a meaningful 
comparison between the power of the different scalars, we 
have chosen in this case higher values for the injected skew- 
ness (S = 0.2) and kurtosis {K = 1.6). The lower part of 
table 121 gives the corresponding mean and dispersion of the 
final skewness and kurtosis obtained from 1000 realistic sim- 
ulations. 



4 RESULTS 

In the present section we have tested the performance of 
the different normalised scalars^ to detect generic devia- 

^ Note that in this section, even if not said explicitely, when talk- 
ing about scalars we will always refer to the normalised quantities. 
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Table 2. Mean and standard deviation of the skewness and kurtosis of maps simulated through the Edgeworth expansion for different 
injected values of S and K . The results have been obtained from 1000 simulations for each considered case. 





injected S 


injected K 


<S> 


o-s 


<K> 


O-A' 




0.00 


0.00 


-1.41 X 10-3 


3.23 X 10-2 


-9.02 X 10-3 


3.38 X 10-2 


Ideal 


0.08 


0.00 


1.18 X 10-2 


3.20 X 10-2 


-8.08 X 10-3 


3.45 X 10-2 




0.00 


0.40 


8.97 X 10-"* 


3.45 X 10-2 


5.82 X 10-3 


3.62 X 10-2 




0.00 


0.00 


2.82 X lO-'' 


4.59 X 10-2 


-1.20 X 10-2 


5.48 X 10-2 


Realistic 


0.20 


0.00 


2.21 X 10^2 


4.63 X 10-2 


-1.17 X 10^2 


5.33 X 10-2 




0.00 


1.60 


-1.55 X 10-3 


4.55 X 10-2 


1.78 X 10-2 


5.21 X 10-2 



tions from non-Gaussianity using the sets of simulations 
previously described. In order to discriminate between the 
Gaussian and non-Gaussian cases, we have applied the 
well-k nown Kolmogorov-Smirnov (KS) test (e.g. lVon Missea 
Il964r for each scalar. Given a scalar map, the test simply 
consists on obtaining the maximum distance between the 
empirically constructed cumulative function of the scalar 
map and the expected theoretical one corresponding to the 
null hypothesis (in our case that the underlying tempera- 
ture map is Gaussian). In most cases this theoretical cumu- 
lative function can be obtained, either analytically or nu- 
merically, for each scalar from the pdf's given in tabled (see 
appendix ^ for the analytical expressions of the cumula- 
tive functions of some scalars). Nonetheless, for two scalars 
(derivative of the gradient and Gaussian curvature) with 
complicated pdf's we found more convenient to obtain it 
simply as the average of the empirically constructed cumu- 
lative function of the 1000 Gaussian simulations. 

For a given normalised scalar, the procedure is as fol- 
lows. First we obtain the scalar map for our sets of Gaus- 
sian and non-Gaussian simulations. Then we construct the 
probability distribution of KS distances for each set of simu- 
lations. Afterwards we compare the Gaussian case with the 
corresponding non-Gaussian one (either ideal or realistic). 
The capability of the scalar to discriminate between both 
cases is quantified by obtaining t he power p o f the test at a 
given significance level a (see e.g. lGowar]ll99?ll . For a signif- 
icance level a, we reject the null hypothesis if a simulation 
has a KS distance, (Iks, higher than that of a fraction 1 — a 
of the Gaussian simulations. Therefore a defines a critical 
value dc below (above) which we accept (reject) the null hy- 
pothesis, p is defined as the fraction of simulations of the 
alternative hypothesis with values of cLks higher than dc- 
Therefore, for a fixed value of a, a large value of p indi- 
cates that there is a small overlap between both probability 
distributions and thus the models can be distinguished. We 
have performed this study for all the considered scalars, as 
well as for the normalised temperature map, using the ideal 
and realistic sets of Gaussian and non- Gaussian simulations 
described in the previous section. The power of the scalars 
to discriminate between the different cases has then been 
compared. 



4.1 Ideal case 

As a first step, we have applied our Gaussian analysis to the 
three types of ideal CMB simulations described in table |2| 
Gaussian, skewness {S — 0.08) and kurtosis {K = 0.4). The 



left panel of figure^shows the probability distribution of KS 
distances, p{dKs), obtained from the temperature maps of 
1000 Gaussian (solid) , skewness (dashed) and kurtosis (dot- 
ted) CMB simulations. As expected from the small values 
of skewness and kurtosis present in the non-Gaussian sim- 
ulations (table I^J, the three curves completely overlap and 
thus the temperature distribution cannot dicriminate be- 
tween the different cases. Similarly, we have performed the 
KS test for each of the normalised scalars. Figure |5| shows 
the KS distance probability distribution of the normalised 
Laplacian, shape index, derivative of the gradient and frac- 
tional anisotropy. As seen in the figure, the Laplacian (top 
left) is able to discriminate very well between the Gaussian 
and both non-Gaussian models, since the probability dis- 
tribution corresponding to the Gaussian case overlaps only 
slightly with the ones of the non-Gaussian simulations. The 
shape index (top right) and the derivative of the gradient 
(bottom left) are able to discriminate between the Gaussian 
and skewness simulations but the power of the test is lower 
for the kurtosis case. For the fractional anisotropy the prob- 
ability distribution of the skewness simulations completely 
overlap with the Gaussian one, whereas the case of the kur- 
tosis simulations is clearly separated. 

Table |5] gives the power of the normalised temperature 
and the different normalised scalars to discriminate between 
the ideal Gaussian and skewness {S — 0.08) simulations for 
two signficance levels (a — 0.01, 0.05). The same information 
is also given for the ideal kurtosis {K = 0.4). We remark that 
the distortion and the shear have the same discriminating 
power since they are related through a strictly monotonous 
function. Thus we present the results only for the distortion. 

As seen in tabled the best scalars at detecting an excess 
of skewness are the fractional isotropy, the Laplacian and the 
shape index, all of them with p > 99 per cent for a — 0.05, 
and also the derivative of the gradient, with a slightly lower 
power. Regarding the kurtosis simulations, the scalars that 
perform better are the Laplacian, the fractional anisotropy, 
the determinant and the curvedness, which all give values 
of the power ^ 99.7 per cent for a = 0.05. Note that the 
performance of the temperature to discriminate between the 
Gaussian and non-Gaussian models is very poor and many of 
the scalars perform better than the temperature, especially 
for the kurtosis case. 

Given the very small amount of skewness and kurtosis 
injected in the temperature distribution of the non-Gaussian 
simulations (see table|2|l, it is not surprising that we can not 
discriminate between the different considered models using 
this quantity. However, we may wonder why some of the 
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Figure 1. Probability distributions of the KS distances p{dxs) corresponding to the CMB normalised temperature maps, each of them 
obtained from 1000 simulations, for the ideal (left panel) and realistic (right panel) cases. The different lines correspond to the Gaussian 
(solid), skewness (dashed) and kurtosis (dotted) sets. 

Table 3. Power, p, of the KS test at two different significance levels {a = 0.05, 0.01) to discriminate between the Gaussian and non- 
Gaussian simulations for the normalised temperature and the normalised scalars. We show the results for both the ideal and realistic 



Ideal case 



Realistic case 



S=0.08 



K=0.40 



S=0.20 



K=1.60 



a = 0.05 a = 0.01 a = 0.05 a = 0.01 a = 0.05 a = 0.01 a = 0.05 a = 0.01 



Temperature 


7.5 


1.9 


8.4 


2.3 


5.3 


1.3 


4.6 


1.5 


Laplacian 


99.8 


98.1 


99.8 


99.0 


99.3 


96.7 


99.9 


98.8 


Greatest eigenvalue 


92.7 


80.1 


92.7 


75.2 


69.7 


41.6 


75.6 


51.2 


Lowest eigenvalue 


92.0 


78.3 


93.1 


76.6 


73.1 


46.9 


78.1 


56.5 


Determinant 


5.9 


0.7 


99.7 


99.2 


4.3 


0.4 


99.5 


98.4 


Distortion 


5.3 


0.8 


73.4 


48.7 


4.3 


0.6 


48.6 


25.7 


Ellipticity 


95.2 


86.9 


61.1 


27.8 


98.4 


91.9 


78.0 


32.9 


Shape index 


99.5 


97.9 


89.3 


63.2 


99.1 


97.1 


96.1 


77.3 


Fractional anisotropy 


6.4 


2.4 


99.7 


98.9 


6.5 


1.8 


100.0 


100.0 


Fractional isotropy 


99.9 


98.8 


91.2 


73.6 


99.7 


97.9 


93.8 


72.4 


Gradient 


8.7 


2.1 


66.7 


43.1 


4.6 


1.2 


72.4 


51.6 


Derivative of the gradient 


98.3 


92.4 


42.5 


10.2 


99.7 


98.7 


74.8 


57.0 


Curvedness 


6.3 


1.2 


99.8 


99.3 


5.5 


0.9 


95.0 


85.6 


Gaussian curvature 


6.8 


1.7 


7.4 


1.9 


7.9 


1.8 


61.2 


37.5 



scalars are able to detect the non-Gaussianity. The reason is 
that this non-Gaussianity can be amplified in certain scalar 
maps. To illustrate this point, we have calculated the mean 
value and dispersion of the skewness of the Laplacian for the 
ideal Gaussian and skewness simulations, finding the values 
(-0.86±4.03) X 10"^, and (2.89±0.41) x 10"^ respectively. 
Comparing these values with those given in the upper part of 
table |5| that were obtained directly from the CMB temper- 
ature simulations, it can be seen that the skewness has been 
significantly enhanced in the Laplacian map. Therefore, for 
the considered case, the non-Gaussianity can be much more 
easily detected using the Laplacian than the temperature 
map. An analogous argument holds for the kurtosis simula- 
tions. In this case we obtained values for the kurtosis of the 



Laplacian map of (0.01±7.87) x 10"^ , and (7.11±0.86) x 10"^ 
for the ideal Gaussian and kurtosis simulations, respectively. 

From the results of table El it also becomes clear that 
one scalar can be very good at discriminating one type of 
non-Gaussianity whereas its performance can be very poor 
for a different type. This points out again the fact that it 
is not a trivial task to find a unique Gaussianity test which 
is optimal in all cases and motivates the investigation of a 
relatively large set of scalars. 

We may wonder how our method compares with 
other Gaussianity tests proposed in the literature. For in- 
stance, techniques based on spherical wavelets have pro- 
vided to be a very us eful tool to study the Gaussian- 
ity o f the CMB (e.g. iBarreiro et al.l 12001 ICavoii et alJ 
I2OOII iMartinez-Gonzalez et alJ I200I fvielva et alJ l2004l) . 
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Figure 2. Probability distribution of the Kolmogorov-Smirnov distances corresponding to the normahsed Laplacian (top left panel), 
the shape index (top right panel), the derivative of the gradient (bottom left panel) and the fractional anisotropy (bottom right panel) of 
the ideal simulations. Each probability distribution is obtained from 1000 simulations for the different considered cases: Gaussian (solid 
line), injected skewness (dashed), and injected kurtosis (dotted). 



In particular, iMartinez-Gonzalez et all i2002l) studied and 
compared the performance of the spherical Mexican hat 
wavelet and the Spherical Haar Wavelet using simulations 
generetated through the Edgeworth expansion. Unfortu- 
nately, that work is not directly comparable to the anal- 
ysis that we present, since there are several differences be- 
tween the used methods and simulations, most notably the 
fact that they use the Fisher dicriminant to separate be- 
tween Gaussian and non- Gaussian distributions and the dif- 
ferent levels of injected skewness and kurtosis. Nonetheless, 
by looking at their results in table 2 we can get some in- 
sight on the comparative performance of the scalars versus 
the spherical wavelets. In particular, it is clear that the best 
of our scalars would outperform the spherical Haar wavelet, 
whereas it seems that they would provide comparable levels 
of detection to those of the spherical Mexican hat wavelet. 



4.2 Realistic case 

We have also applied our technique to the realistic CMB sim- 
ulations described in section |3] As in the ideal case, we have 
considered three sets of simulations (see table I^J: Gaussian, 
injected skewness (5* = 0.2) and injected kurtosis [K — 1.6). 
As already explained, we have increased the levels of non- 
Gaussianity in these simulations to allow for a better com- 
parison between the different scalars in the presence of noise 
and a mask. 

As for the ideal case, we find that the probability dis- 
tributions of KS distances, p{dKs), corresponding to the 
Gaussian and non-Gaussian simulations completely overlap 
for the normalised temperature (see right panel of Fig. 0. 
Regarding the scalars, we give in Fig. |3 the KS distance 
probability distributions of the normalised Laplacian, shape 
index, derivative of the gradient and fractional anisotropy, 
for the three sets of realistic simulations. The results are 
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1000 - 
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Figure 3. Probability distributions of the KS distances, corresponding to the Laplacian (top left panel), the shape index (top right 
panel), the derivative of the gradient (bottom left panel) and the fractional anisotropy (bottom right panel) of the realistic simulations. 
Each probability distribution is obtained from 1000 simulations and the different lines correspond to the Gaussian case (solid line), 
injected skewness (dashed) and injected kurtosis (dotted). 



qualitatively similar to those of the ideal case: the Laplacian 
can distinguish reasonably well the Gaussian case from both 
types of non-Gaussian simulations, the shape index and the 
derivative of the gradient are good at detecting the skewness 
case whereas the fractional anisotropy can discriminate only 
between the Gaussian and kurtosis simulations. 

The results regarding the power of the test at two dif- 
ferent levels of significance (a = 0.01, 0.05) for the tempera- 
ture and the considered scalars are given in tablej^for both 
the skewness and kurtosis realistic cases. The comparative 
performance of the different scalars is very similar to that 
obtained for the ideal case, although some scalars seem to 
be more affected by the noise than others. In particular, 
for the skewness case, the fractional isotropy, the Laplacian 
and the shape index give again high values of the power, but 
they are outperformed by the derivative of the gradient, that 
seems more robust under realistic conditions. Regarding the 
kurtosis simulations, the fractional anisotropy, the Lapla- 



cian and the determinant perform again very well, whereas 
the curvedness is more affected by the presence of noise and 
gives a somewhat lower discriminating power. 



5 CONCLUSIONS 

M05 introduced a set of normalised scalars constructed on 
the sphere from the covariant derivatives of the tempera- 
ture field. The possible use of these quantities for a Gaus- 
sian analysis of the CMB was also proposed. In the present 
work, following M05, we have presented a Gaussian analysis 
for CMB data based on the statitiscal properties of the nor- 
malised scalars. The comparative performance of the differ- 
ent scalars has been tested using Gaussian and non-Gaussian 
simulations generated through the Edgeworth expansion. In 
particular, we have used two different types of non-Gaussian 
simulations: with injected skewness and with injected kur- 
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tosis. In addition we have considered an ideal case (where 
noiseless all-sky simulations were constructed) and a realis- 
tic case (containing anisotropic noise and a mask). 

In order to quantify the power of the different scalars 
to detect non-Gaussianity we have used the Kolmogorov- 
Smirnov test. We find that, in most cases, the scalars am- 
plify the non-Gaussianity present in the temperature map 
and produce higher detections than those obtained directly 
with the temperature. In particular, the best scalars to dis- 
criminate between the Gaussian and injected skewness cases 
are the Laplacian, the fractional isotropy, the shape index 
and the derivative of the gradient. For the kurtosis case, 
the highest powers are found for the determinant, the frac- 
tional anisotropy, the Laplacian and the curvedness. Note 
that some scalars can be very good at detecting one type of 
non-Gaussianity but perform very poorly for discriminating 
the other type. 

In future works, we expect to test the discrimi- 
nating power of these quantities using physically moti- 
vated non-Gaussian models such as the quadr atic potencial 
model IILiguori. Matarrese fc Moscardinal2003i) and also ap- 
ply them to the analysis of the WMAP data. 
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APPENDIX A: CURVEDNESS 

In this appendix we describe a new scalar, the curvedness 
l|Koenderinkl Il99n i . that was not included in the study of 
M05. The curvedness, c, is a scalar defined in terms of the 
eigenvalues of the negative Hessian matrix of the field Ai 
and A2 by the following expression: 



lyXfTAf 



(Al) 



By construction c is a positive quantity and is closely related 
to the curvature of the field. Note that points with larger 
values of Ai or A2 (positive or negative) will have larger 
values of c whereas it becomes zero for fiat areas in the 
initial field (i.e. Ai = A2 = 0). 

Similarly to the other scalars, the normalised curved- 
ness is constructed replacing Ai and A2 by the normalised 
eigenvalues in equation <AH : 



^=2V^i+^^ 



(A2) 



When the initial field is Gaussian distributed, the probabil- 
ity distribution function of the curvedness is given by: 






(A3) 



which is independent of the power spectrum of the field. D 
is the Dawson's function, defined by the integral equation: 



D(x) 



e dt 



(A4) 



The correlations of the normalised curvedness with the other 
normalised scalars are given in table IXT1 



APPENDIX B: FRACTIONAL ANISOTROPY 

In this appendix we descr ibe another interesting scalar, 
the fractional anisotropy fa JBasser fc Pierpaol i"l996^ . This 
quantity has been used in the analysis of medical images, 
including the inv estigation of some n eural diseases like as- 
trocytic tumors llBeppu et alJ l2003r or multiple sclerosis 
JWerring et al.lll999h . 

fa is defined as a function of the eigenvalues Ai and A2: 



fa = 



Al — A2 



V2,/XiTX: 



(Bl) 



which is proportional to the ratio of the distortion over 
the curvedness. By construction fa is dimensionless and 
bounded, < /a < 1. It encodes information about the 
asymmetry between the eigenvalues of the field. 

The normalised fractional anisotropy fa is given by 



fa 



1 Al - A2 

v2^Afn 



(B2) 



When the initial field is Gaussian, fa will follow a probability 
density function given by 



Pifa) = 



2fa 



1-/2 ni + /I 



(B3) 



which, by construction, is independent of the initial power 
spectrum of the field. 

The correlations of fa with the other normalised scalars 
are shown in table IXTl 



APPENDIX C: FRACTIONAL ISOTROPY 

The fractional isotropy fi is proportional to the ratio of the 
Laplacian over the curvedness: 

1 Al -f A2 



/» = 



V2,/JfTX: 



(CI) 



By construction fi is dimensionless and bounded, — 1 < /i < 
1. 

Following the usual construction, the corresponding 
normalised scalar fi is given by: 



/. 



Ai-f A2 



For an original Gaussian field fi follows the pdf: 

p(/0 



i^-fiy 



(C2) 



(C3) 



which is independent of the power spectrum of the field. 

We include in table lAll the correlations of the nor- 
malised fractional isotropy with the rest of normalised 
scalars. 



APPENDIX D: CUMULATIVE FUNCTIONS OF 
THE NORMALISED SCALARS 

For some scalars, it is possible to obtain an analytical expres- 
sion of its corresponding cumulative function by integrating 
the pdf's given in table which are valid for the case of an 
initial Gaussian temperature field. 

In particular, in table IHTl we give the cumulative func- 
tions of the normalised Laplacian, shear, distortion, ellip- 
ticity, shape index, gradient, fractional anisotropy and frac- 
tional isotropy. 
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Ai A2 



A4 



y 



X- 



Da 



i^G 



fa 



U 



c 0.26 -0.26 
fa 0.23 -0.23 
fi 0.83 0.83 0.92 



0.26 0.61 0.62 

-0.77 0.45 0.54 

-0.99 



0.21 1.00 -0.25 

-0.69 -0.25 1.00 

-0.53 1.00 



Table Al. Correlations of c, fa and fi with the rest of normalised scalars (obtained from Gaussian simulations). 



Normalised scalar 


Notation 


Domain 


cumulative 


Laplacian 


A+ 


(-00,00) 


F(A+) = i(l + e./(^)) 


shear 


y 


(0,00) 


F(y) = 1 - e-4i/ 


distortion 


X- 


(0,00) 


F(A_) = l-e-^- 



cUipticity 



(—00, 00) -P(e) 



1 1 

2 y/T+SS 



2 v^T+Sp" 



e < 
e > 



shape index 



(-1,1) F{1) 



1- 



^l + 2cot2(fj) 

1 

^l + 2cot2(fj) 



t< 

;;> 



gradient 


9 


(0,00) 


F{g) = l-e-~a 




fa 


(0,1) 




fractional anisotropy 


nW-1-/^ 


fractional isotropy 


h 


(-1,1) 


PI f ^ 1(11 fi 





Table Dl. Cumulative functions of some normalised scalars, for the case of an initial Gaussian field. 



